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Abstract 

Angular perturbations modify the band structure of armchair (and other metallic) carbon nan- 
otubes by breaking the tube symmetry and may induce a metal-semiconductor transition when 
certain selection rules are satisfied. The symmetry requirements apply for both the nanotube 
and the perturbation potential, as studied within a nonorthogonal vr-orbital tight-binding method. 
Perturbations of two categories are considered: an on-site electrostatic potential and a lattice defor- 
mation which changes the off-site hopping integrals. Armchair nanotubes are proved to be robust 
against the metal-semiconductor transition in second-order perturbation theory due to their high 
symmetry, but can develop a nonzero gap by extending the perturbation series to higher orders or 
by combining potentials of different types. An assumption of orthogonality between vr orbitals is 
shown to lead to an accidental electron-hole symmetry and extra selection rules that are weakly 
broken in the nonorthogonal theory. These results are further generalized to metallic nanotubes of 
arbitrary chirality. 

PACS numbers: 73.63.Fg, 73.22.-f, 61.46. +w, 71.15.-m 
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I. INTRODUCTION 



The subject of metal-insulator transitions has been studied for decades.— It is well under- 
stood that a metal-insulator transition is typically related to the breaking of a specific sym- 
metry of the system, and carbon nanotubes are particularly interesting to study, due to their 
low dimensionality and special helical symmetry. In this paper we investigate the symmetry 
breaking in single-wall nanotubes (SWNTs) due to an external potential which depends on 
the angular coordinate along the SWNT circumference. This perturbation may induce a 
transition in a SWNT changing the type of its electronic structure by opening/ closing the 
band gap in the metallic/semiconducting nanotube .^1^1^1^1^1^ It is important for applications 
in which such gap engineering can potentially be controlled locally, for instance by the field 
of a sharp tip, by a small molecule or by a local gate. 

SWNT lattice symmetry depends on two parameters, diameter and chirality, which deter- 
mine the type of band structure.^^ About one-third of possible SWNTs are metallic, with one 
dimensional energy subbands crossing at the Fermi level, as confirmed by experiments.— In 
this study, we consider metallic nanotubes and focus mainly on the special case of armchair 
SWNT (A-SWNT), which has a higher symmetry prohibiting an energy gap at the Fermi 
levelii for typical non-chiral perturbations such as stretching, uniform electric field, impu- 
rity potentials and many-body interactions. The same perturbations would open a small 
"secondary" band gap in metallic nanotubes of different symmetry. The high symmetry of 
A-SWNTs is also responsible for the absence of back scattering in the conduction channels.— 
Such ballistic transport is very attractive for future electronic applications^iiS and a method 
to control the conductance of A-SWNTs would be particularly desirable. Different perturba- 
tions have been attempted to modify the electronic structure of A-SWNTs . ^i^i^i^i^i^i^i-^^i-^^i-^^i-^^ 
Our goal is to demonstrate, using symmetry arguments, whether a particular perturbation 
can open a gap at all and how the gap depends on the magnitude of the perturbation 
potential. We will show below that, with minor exceptions, this cannot be a linear depen- 
dence. We call this transition a "metal-semiconductor transition" (MST) since the gap size 
is smaller than in typical insulators.— 

A nonorthogonal tight binding (TB) approach is used to model the SWNT electronic 
structure. Despite its simplicity, the TB approach may include as much important physics 
as more sophisticated models with the right choice of empirical parameters.— In addition. 
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it possesses a great advantage for analytical derivations. Combining the TB approach and 
the summation of perturbation series with a group theory technique, it was shown in previ- 
ous worki^ that mirror symmetry breaking is a necessary condition to mix the two crossing 
subbands (vr and tt*) and open a band gap in A-SWNTs. Here, we find several additional 
significant results: (1) Due to the high lattice symmetry, the second order contributions 
always cancel out and no second order gap opening occurs. One notes that the first or- 
der process is suppressed unless very specific selection rules are satisfied (see below). (2) 
Under potentials of a single angular Fourier component, V^cosg^, the lowest contributing 
coupling order between vr and vr* bands, /iq, is determined by the angular momentum of the 
potential, q, and the index n of the (n, n) A-SWNT as /io = 2n/gcd{2n,q), in which gcd 
is the greatest common divisor. The band gap opening is proportional to |Vq|'^° for small 
perturbation |Vg| <^ vp/R, where vp and R are the nanotube Fermi velocity and radius. 
In a typical experiment, when the perturbation has a small angular momentum (g ~ 1), 
the coupling order is high {^q ~ n) and the gap is small if any. To observe a linear effect 
(/io = 1), a high-g potential must be applied (g = 2n). The high coupling order can be 
reduced by choosing combinations of several angular Fourier components or different types 
of perturbation. (3) Additional symmetry of a particular model may lead to extra selection 
rules for the band gap opening. For example, gaping is forbidden for A-SWNTs with even 
n if perfect electron-hole symmetry is assumed, as in an orthogonal basis. When the slight 
asymmetry between the conduction and valence bands is included, a gap proportional to 
the asymmetry parameter occurs. (4) Significant changes in the A-SWNT density of states 
(DOS) is observed even when the gap is absent. Modification of the low energy band struc- 
ture is mostly determined by the second order perturbation. The DOS is enhanced near the 
Fermi level and simultaneously vp decreases. Peaks of the first pair of van Hove singularities 
are brought closer, resulting in a smaller excitation energy between these subbands. 

Our study is restricted to the case of potentials which are uniform along the tube axis, 
but certain results are easily generalized for perturbations with even/odd axial dependence. 

This paper is organized as follows. We first formulate the model and introduce the 
interaction matrix elements between TB wave functions in Sec. |H] Using nearly degenerate 
perturbation theory (Appendix p) and symmetry-based selection rules, we derive analytically 
the coupling between vr and n* subbands of A-SWNTs for both scalar potentials (Sec. IIII|) 
and tensor potentials (Sec. lIVjl . Comparisons are made with numerical results from TB band 
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structure calculation. These results are further extended to metallic nanotubes of arbitrary 
chirality in Sec. IVII Finally, we summarize our results in Sec. I VI II 



II. MODEL FORMULATION 



A. Perturbation series in the TB description 



The electronic states are obtained within a nonorthogonal single vr-orbital TB method 
by solving the stationary Schroedinger equation 

Hm, = E,S^,, (1) 

where H and S are the Hamiltonian matrix and overlap matrix, respectively. For the sake 
of simplicity, only nearest neighbor hopping integral 70 = —3.033 eV and overlap integral 
= 0.129 are considered.— 

The wave function of an unperturbed A-SWNT can be expressed as a linear combination 
of the two periodic functions M§(k) = '''''*~'"V(j' ~ J'ic)' where ^ = A^B label the 

two sublattices and (p{Y — Vi^) is the atomic orbital function localized at rj^: 

A — 1 

where a = ±1 denote the conduction and valence bands, ta's are the nearest neighbor bond 
vectors and we refer to A = 1 as the circumferential direction of A-SWNTs in the following. 
The wave vector k is composed of an axial component kt and a quantized circumferential 
component = m/R, with m the angular momentum. 2(j)„(k) = arg[— cr/(k)] indicates the 
phase difference of the coefficients before M^(k) and ^^(k), which make a pseudo-spinor (see 
also Sec. lVI|) . 0o-(k) is constant through the whole range of kt for the two crossing subbands: 
4>n{kt) = it/3 and 4>TT*{kt) = (p-K — 7r/2. Clearly, the n and vr* subbands are orthogonal due 
to this phase difference. 

Consider a perturbation Hi which is uniform in the axial direction, then kt is conserved. 
To obtain the low energy behavior of vr and vr* subbands, we write the effective 2x2 per- 
turbation Hamiltonian matrix using nearly degenerate perturbation theory (see Appendix): 
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where the matrix element Hapikt), with a{(3) = tt or vr*, can be represented by the sum of 
the perturbation series over all possible coupling orders fi as 

n (-^.) 

One notes that all intermediate states ii = 1 ... /i — 1) are different from \E'^ and "^.j^* 
by definition. Figure ^ illustrates an example of second order coupling between vr and vr* 
subbands via four different paths. Also shown is the phase angle of the intermediate states, 
(f)±{kt,n ± q), relative to (p-,, and (p-,,*- The relation between these phase angles will be 
discussed in detail in the following sections. 

B. Interaction matrix element between Bloch states 

Below, we derive the interaction matrix element (^| Hi |\t'') of the angular perturbation 
of a single Fourier component. Hi = \q cos q{6 — Oq), where Vg could be either scalar or 
tensor, is defined as the minimum angular displacement between the vertical mirror 
planes (or glide planes) of the A-SWNT and the mirror planes of the potential.— Assuming 
that {(fi \ Hi \(f2) is nonzero only when (fi and (f2 are centered on the same atom or nearest 
neighbor atoms, the interaction matrix element can be decomposed into an on-site term and 
an overlap term: 

(^i)on-sitc = (v^(r)|V,e'«^|vp(r))cos(0-0') 

(^i)overiap = E (^(^ " ^^2)1 V,e^<^^ |^(r + v,/2)) COS (0 + 0' - ^^^^A " hz!) (5) 

A=l ^ ^ 

where it is used that y9(r) is real and invariant under 9 ^ —9 inversion. 5m-m.',±q arises from 
the conservation of the angular momentum. This conservation law also allows q to differ by 
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-f^TTTT* {kt 



(3) 



multiple of 2n, an angular analog of the reciprocal lattice vector, which is not considered 
here (see discussion in Sec. IVI|) . The interaction matrix can be further simplified depending 
on the type of Vg: 

Scalar Potential: =const. Define {^p{r)\ V^e**^^ lv^(i")) = Uq, then 

{^{r-r,/2)\Vqe^^<^Mr + r,/2)) ^ ^0 [(^(r + r,/2)| V^.e^"^ |^(r + r,/2)) 



+ (y.(r-r,/2)|Ke^''^|^(r-r,/2))] 
2 



0Uq COS 



(6) 



(^l)ovcrlap = I^Uq 



A=l 



im'Bx+iktZx ^ 



+ C.C. 



-^0f/, COS (0-00 {a\f\+a'\f'\), 



and Eq. is reduced to 

(-^l) on-site ~l~ (^l) overlap 



^ Uq C0S( 



l--0(a|/|+a'|/' 



(7) 



(8) 



Tensor Potential: The Fourier component of a tensor potential can be expressed in the 
second quantization formalism as 

A=l (»j)A 

where pairs are confined to first nearest neighbors with bonds along direction 

and S'jx^q is the corresponding change of the hopping integral. The on-site term is 
absent while the overlap term is reduced to: 

3 

(^i)overiap ~ ^A(/i;t; m', a) 



A=l 
3 



E Sjx,g COS ( <; 
A=l ^ 



m + m' 



Ox - hzx . 



(10) 



Comparing Eq. © and Eq. (jlUp . one can conclude that the interaction matrix elements from 
a scalar potential and a tensor potential have quite different dependence on the phase angle 
(f) and the quantum numbers kt, m. We will show in next section that certain selection rules 
can be derived for scalar potentials of general form and also for simple tensor potentials 
allowing summation over A. 



III. SCALAR POTENTIAL 



Assume that a scalar perturbation in the form of Hi = Vq cos q{6 — 60) is apphed to the 
{n, n) A-SWNT. Using the interaction matrix element from Sec. IHl we can now derive the 
/i-th order perturbation matrix elements within nearly degenerate perturbation theory: 



AfiqOo 



Qim) 



JJcos( 



J=l 



n 

i=l 



i-1 



n i-E?) 

i=l 

cos(0^_i + Afiqn/Sn 



(11) 



where subscripts "0" and "/x" correspond to the initial state \1/q, and final state \l//3 respec- 
tively, with mo = = n. We stress that Ef = E^.{kt,mi) = —ai'-yo\f{kt,mi)\, since 
the factor (1 — 0o'i\fi\)~^ is canceled by those from the wave functions. Pais is the total 
phase of the perturbation term and corresponds to inner products of the pseudo-spinors 
(see Sec. IVII) . Q corrects for contributions from a nonzero nearest neighbor overlap. The 
intermediate states = \E'o--(A;(, mj)} satisfy the conservation law of angular momentum 
and have the following constraints: 



mi_i -nii = ±q, i = 1, . . . ,fi - 1 
m^_i — = ±q + multiple of 2n, 



(12) 



which leads to ^f^^ ('^i-i 



nii 



A/ig + multiple of 2n = 0. The role of a nonzero A/i can 



be seen from the extra phase factor e~*^^''^° in and the corresponding term Afiqir/Sn in 
Pais- Direct evaluation of Eq. (Illj) with all possible sets of {^^j} is formidable, nevertheless, 
one can get useful information by applying symmetry arguments. 



A. Off-diagonal coupling: gapping of A-SWNTs 

We first study the off-diagonal term -ff^^l and replace with irii = 2n — m in Eq. (Illj). 
which is allowed by conservation of angular momentum. The energy denominators and the 
function Q remain unchanged because \f{kt,m)\ = \f{kt,m)\, while the sign before Afiq 
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changes. Notice that for A-SWNTs, f{kt,m) = e''^''/^f*{kt, m), so that by the definition of 
0, cos(0 — (f)') = aa' cos(0 — (p'). Define = ^„i{kt, 2n — rrii), then 

= a^a^,P^^,{{^,}) = -P^^.({^,}), (13) 

which leads to 

HirMm) + Hi'^Mm) ^ (^^J sin{A^,qe,)P^Am)- (14) 

Since the above relation is true for all possible sets of intermediate states at all coupling 
orders, one can conclude that: (a) The coupling between vr and vr* subbands is always zero 
if 6q = 0, i.e., when a mirror plane of the potential overlaps with one of the vertical mirror 
planes (or glide planes) of the A-SWNT. This is an explicit result of the mirror symmetry 
requirement.^ (b) The coupling is zero if A/i = 0, which results from refiection symmetry of 
the energy bands: Ecr{kt, m) = Ecr{kt, 2n — m). This excludes the possibility of any nonzero 
second order contribution, i.e., fi = 2, AyU = 0. In other words, a second order band gap is 
forbidden in A-SWNTs. (c) The next lowest possible A/i satisfying the angular momentum 
conservation in Eq. (fT^ is fiQ = 2n/gcd{2n,q), which is also the lowest contributing order 
of the perturbation series.— For small Uq, the band gap opening will be at least the fiQ-th 
order in Uq, 

Eg ^ 2\H^^,{kF) \ ~ ^u^" sm{fioqeo)h{q,n) + (terms of /i > fio), (15) 
K 

where kp = 27r/3a is the Fermi point with a ^ 2.5 A. u = UqR/v-p is the dimensionless 
potential and h{q,n) is a complicated function depending on the angular momentum of the 
potential, q, as well as the A-SWNT index, n. 

The summation over all possible intermediate states can be further simplified by com- 
bining the original process with {\E'i} and the reversal process with {\E'f = '$_^^_.{kt,2n — 
m^„j)}. The sign change of the energy results in an extra factor of (—1)^"^ in the denomina- 
tor and the function Q changes accordingly. It can be proved that P^^*({\Eff }) = P^^. ({\['j}) 
and therefore at small the relation holds: 

i=l 
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Within the orthogonal model (0 = 0), only the first term in Eq. ()16p exists and is nonzero 
when /io, and therefore /i, is odd. This constrain on fio results from the invariance of the 
inner product of pseudo-spinors upon reversal operation, in combination with the electron- 
hole symmetry E^_^{kt,m) = —E'^{kt,m). The latter, however, is not an intrinsic property 
of SWNTs, but rather due to the nearest neighbor approximation. For instance, the energy 
band symmetry is broken when the second nearest neighbor hopping integral is included, 
or, equivalently, if 7^ 0. In the latter case, Q({\E'i}) and — (5({\l'f }) don't cancel out and a 
non-zero band gap proportional to opens for even /iq. Figure |21 plots the variation of band 
gap at g = 2 calculated within the orthogonal and nonorthogonal TB models. At s = 0, the 
(6,6) A-SWNT remains metallic, because the coupling order /xq = 6 is forbidden. At nonzero 
0, a small band gap occurs and increases as a power law of u. It is also found that the band 
gap grows linearly with the magnitude of 0, consistent with the prediction of Eq. (fTBj) . In 
contrast, the band gap curve of a (5,5) A-SWNT only shows a slight increase at nonzero 0, 
because the corrections is of the order of 0^. 

B. Diagonal coupling: renormalization of the Fermi velocity 

Except for a few special cases, for instance with q = multiples of 2n, the coupling order 
between vr and tt* subbands is about the same order of n and the resulting band gap remains 
small. For example, /iq = 2n at g = 1, and fio = n at q = 2. However, the diagonal coupling 
matrix elements Ht^^^ and 11^^*^^* are not necessarily small. Same symmetry arguments can 
be applied here. Upon reflection or reversal operation, P^rTr (or Pn*n*) remains the same and 
one obtains 

Hm^^}) + Hi^2im) cc (^^J cos{Auqe,)PUm) (17) 

OC [p„„ - {-lyppp] - [Pa. + {-lypM ^^i/^i + o{0'), (18) 

1=1 

where a{(3) = tt, tt* and a ^ (3. Since cos(Az/g0o) is always unity whenever Az/ = , the 
lowest contributing coupling order is therefore v = 2. So unlike the off-diagonal coupling, 
nonzero diagonal terms H-j^t^ and Ht^*-,^* can always be obtained from low order angular 
perturbation. In an orthogonal basis (0 = 0), only the first term in Eq. ()18|) remains, which 
corresponds to an energy shift for n and tt* subbands in the same direction when z/ = odd 
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and the opposite direction when v = even. If 7^ 0, relative shift between the two subbands 
always occurs. The values of H^^t, and do not contribute to the band gap opening of 

the A-SWNT, but may influence the Fermi point position as well as the DOS near the Fermi 
level. Assume that = 0, then the second order perturbation summation is reduced to 

'Uq\^ cos^ [0^ - (f)+{kt,n + q)] - cos^ [(f)^* - (f)+{kt,n + q)] 



V H'-^^ = 2 r — ^ cos^ [0^ - (f)+{kt, n + q)]- cos^ [0, 
£1 V2y -E+{kt,n + q) 



' {^^} 

7^/, N 1 cos(2q'7r/3n) ,, , , , , ^ / . , 9\ 

^ ^ ~ -^hT— ^-TT^ + o • 2 /' / - a + O Afc2 . 19 

2 [1 + 2 cos(g7r/n)J 8sm (g7r/2n) ^ ^ 

The TT subband is a decreasing function of \kt\, and since F{kt,q) is an increasing function 
of \kt\ at small g values, it becomes flattened near ±kp as a result of the second order 
perturbation. Similar trend can be found for the ir* subband. The new Fermi points move 
toward kt = as F{±kF, q) < 0, and the renormalized Fermi velocity is given by 

where the prefactor (l+u^/2g^)~^ is due to the normalization of the perturbed wave function. 
vp can also be estimated from a chiral gauge transformation^^ as Jo{2u / q)vF ~ (1 — u^/g^)fF, 
with Jo the Bessel function of the first kind. The renormalized Fermi velocity of a (10, 10) 
A-SWNT is plotted in Fig. IHl and excellent agreement is found between the TB results and 
the analytical predictions. Under potentials of small g's, the shape of vr and vr* subbands is 
strongly perturbed and the low energy DOS is greatly enhanced (see insets of Fig. El), which 
becomes more evident for large radius A-SWNTs even at a relatively weak perturbation. 

One may note that, for A-SWNTs, the finite curvature shifts the Fermi point further 
toward k = while inclusion of a-orbitals may also modify the magnitude of the band gap 
opening . ^1^1^^ For the radius range considered in this paper, the correction remains small, 
and since our symmetry arguments are based on the lattice geometry of A-SWNTs (mirror 
or angular symmetry), the selections rules are not affected. 
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IV. TENSOR POTENTIAL: q = 2 



Another type of perturbation is realized by changing the hopping integrals between neigh- 
boring atomic orbitals, for example, by elastic deformations. Due to the high flexibility, 
carbon nanotubes can sustain remarkable deformations, which cause drastic changes in the 
electronic properties. Stretching, twisting and squashing the nanotubes have been attempted 
both theoreticall y^i^i^i^^i-^^i-^^i^^i^^ and experimentally.— Here, we apply a radial defor- 
mation across the A-SWNT and derive the condition to open a secondary band gap. 

Assume that a uniform stress is applied along y direction, which is rotated from the A- 
SWNT's vertical mirror planes (or glide planes) by ^o- The cross section of the A-SWNT is 
distorted into an ellipse with the two axes given by 

Ry = RQ{l-e), R, = Ro{l + ue), (21) 

where e is the strain along y direction and u is the cross section Poisson ratio. At small 
strain, the hopping integral between nearest neighbors can be linearly expanded as 

- 70 ~ (^^^ Arx,i = S-fx (^^^ + cos2^A,i^ (22) 

,,,._2,.„,1,.,(^.4)^(|I)^^. (23, 

where 6\^i corresponds to the midpoint of the i-th bond along r;^ direction. The first term in 
bracket in Eq. (j22j) has no angular dependence and does not break mirror symmetry because 
6*2 = ^3 and thus ^72 = 573. Since the Poisson ratio was found to be close to unity,— the 
second term dominates, which corresponds to an angular momentum q = 2. 

The perturbation matrix element Haf3 can be expanded into a perturbation series in a 
similar fashion as in Eq. (fTTj). but now with Qa/sii'^i}) = 1 and 

M 3 

■,mi_i,ai-i;mi,ai) (24) 

i=l A=l 

where g\ is defined in Eq. dTUI) . Using the fact that ^72 = ^73, one can prove that 

3 3 
'^gx{kt;m,a.m',a') = aa'^gx{kt;2n - m,a.2n - m ,a'), (25) 

A=l A=l 

3 

— ''^^gx{kt;2n — m,—a';2n — m,—a). (26) 



A=l 
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By applying the reflection operation ^ on the intermediate states, Eq. (fTSjl and ()14p 
are recovered, as required by mirror symmetry conservation. The lowest contributing order 
is therefore /io = 2n/gcd(2n, 2) = n. Now apply a reversal operation ^> — >■ as in Sec. IIIII 
Since E^^{kt, 2n — m) = —E^{kt, m), from Eq. ()2(j|l one has 



(-11^ 



0. (27) 



This means that a band gap does not occur in any perturbation order and such strain 
cannot induce MST in A-SWNTs. Certainly, a hidden symmetry, namely, the electron-hole 
symmetry E^^{kt, —m) = —E^{kt,m), forbids the vr and ir* subband mixing. By including 
the second nearest neighbor interactions, this symmetry can be weakly broken and a finite 
band gap occurs.— The magnitude of this band gap depends strongly on the parity of the 
A-SWNT index n. It was reported earlier that squashing a (6, 6) or (8, 8) A-SWNT does not 
induce a MST in the range of elastic deformation , and the nanotube remains metallic 
until the two opposite walls are brought close enough to form new bonds. We prove that the 
vanishing band gap is caused by the high coupling order between vr and vr* subbands, fiQ = n, 
and the additional smallness of the overlap integral and higher neighbor interactions, which 
make it impossible to observe the MST effect until the A-SWNT collapses. The situation 
is quite different, for example, for a (5, 5) A-SWNT which has an odd number index and a 
smaller coupling order, and a finite band gap was observed at moderate deformation.^^ On 
the other hand, when the radial deformation is large enough to induce strong tt-ct interaction, 
the single 7r-orbital description may no longer be sufficient. 



V. COMBINATION OF DIFFERENT TYPES OF POTENTIALS 

One way to reduce the coupling order fiQ is by combining potentials of different angular 
momentum. For example, we have shown that by applying a scalar potential of the form 
of Vo (sin 6' + sin 2^), /io can be reduced to three for all values of n.— The combination of 
elastic radial deformation and uniaxial electrostatic potential will have a similar effect. By 
choosing appropriate angular momentum and relative position of the two components, the 
coupling order can be even lowered to fiQ = 2, as shown below. 

Assume that a scalar potential (denoted as U) is applied on an A-SWNT together with 
a tensor perturbation (denoted as £^ oc 6'~fq) of same angular momentum q, but with an 
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angular difference 9d between the mirror planes of these two components. As shown in 
Sec. mil and Sec. IIVI the second order contribution from either component is zero, but the 



cross terms do not necessarily vanish. For the second order coupling between vr and vr* 

£ U US 

subbands, — ^ ^^(fct, fn) — ^ there are 8 different cross terms with m = n ± q and 
0" = ±1. At kf = k-p, the cross terms adds up to 

3 

E H^XrossikF) = sin{qe,) YjoTZ^'Tl^l [20+(fcp, n + q) - {n + q/2)e,] 



- i\El{kp,n + q)\ 
I To I 

oc uesm(q6d) , q <^ n, (28) 

q 

where u = UqR/vp and e ~ ^7/70 are the dimensionless potential and strain respectively. 
According to Eq. ()28|1 . the coupling between vr and vr* subbands is largest when 6d = vr/2g 
and always vanishes whenever the mirror planes of U and S overlap. The magnitude of the 
band gap, Eg ^ 2| ^i^^^l ctossI^f)!, is linear in both u and e, and the dependence on the 
A-SWNT radius R (or index n) is very weak. This differs from the situation of mixed scalar 
potentials of different angular momentum, in which case Eg decreases with R as an inverse 
power law.— 

In the case of radial deformation (Sec. IIV|1 . the tensor perturbation on the hopping 
integrals and the effective on-site potential^ have overlapping mirror planes so that no 
second order coupling occurs. Assume that one applies on the A-SWNT a scalar potential 
with the q = 2 component shifted 6^ = 7r/4 relative to the stress, for example by changing 
the electrostatic environment around the A-SWNT. The resulting band gap is plotted in 
Fig. E] for A-SWNTs of different radius as a function of u. The hopping integral under 
deformation is assumed to change as 7 oc r~^, where the new bond length r is calculated 
from Eq. ()21|) . The Poisson ratio is taken to be unity. The numerical values of Eg clearly 
follow a linear dependence on u and e when the perturbation is weak. The radius dependence 
is barely seen even at large u, which means that one can always generate a substantial band 
gap in a large radius A-SWNT using only a moderate external potential. Recently, it was 
found that the conductance of a carbon nanotube can be controlled by tuning the voltage of 
the local gate placed near mechanical defects on the tube (kinks or bends). This scenario 
is reminiscent of the combination of scalar and tensor potentials of g = 1, and might be 
related to the resulting second-order band gap. 
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VI. METALLIC SWNTS WITH ARBITRARY CHIRALITY 



In this section, we generalize our derivation of subband coupling to arbitrary metallic 
SWNTs by expanding the TB wave functions near the Fermi point. Only angular pertur- 
bations will be considered here, i.e., the axial wave vector is always conserved. An analogy 
is made to the cancellation rules of the back scattering process in SWNTs.-i For the sake 
of clarity of the derivation, the overlap integral is assumed zero. We remind that some 
cancellation rules will be weakly broken if the electron-hole symmetry is lifted, for example, 
by 7^ 0. 

First, the electronic states are approximated as the product of a plane wave part and a 
pseudo-spinor part by expanding the wave vector near the Fermi Point K of two-dimensional 
graphite: 

^.(k) = ^ I J , 20.(k) = -a^ + arg (4 + %k) + r/, (29) 

where k = (fc^, k^. = rh/R) is measured from K and rj is the chiral angle. It can be proved 
that the definition of ^.(k) here is consistent with that of 0o-(k) in Sec. |n] In contrast to 
Ref. 0, the phase difference due to the sign of a is absorbed in the definition of ^.(k) so that 
the product of two pseudo-spinors is always real. In addition, the wave functions of the two 
crossing subbands of metallic SWNTs remain continuous as a function of kf when the wave 
vector passes through the Fermi point. 

Assume applying an arbitrary angular scalar potential V{6) = VqC'^'^^, where Vq is the 
angular Fourier component of the potential: Vq = (27r)^^ ^ d6V{6)e~^'^^ . The direct coupling 
matrix element between two states ^^^(k) and \E'o-2(k'), with kt = k[, can be reduced to 

M(k, a; k', d') = {k,(T\V{e)\k',a') 

~ ^ ^ Urh—rh'+Q 

Q 

X ^^Q(eoA+0oB)/2 [^^^' + Q^e^B - eoA)/2] , (30) 

where 6oa and Oqb are angular coordinates of any A and B atoms. Angular quantum num- 
ber Q corresponds to the angular analog of the reciprocal lattice vector of two-dimensional 
graphite and accounts for contributions from the short wave-length component of the po- 
tential. When all nonzero Q's are neglected, the matrix element M(k, a; k', a') is reduced to 
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the product of the Fourier transform of the potential and the inner product of two pseudo- 
spinors, comparable with the on-site term in Eq. (0) for A-SWNTs. 

If the perturbation is a tensor potential, i.e., affecting the off-site hopping integrals instead 
of the on-site energies, one can similarly write the coupling matrix as: 

M(k, a; k', a') = {k,a\6^°P \k' 



A=l 



X COS 



+ 



k + k' 



K TA 



(31) 



where S'jx^g is the discrete angular Fourier transform of the change in hopping integrals: 
^1\,q = Comparing with Eq. one finds that matrix elements for 

the two types of perturbations have different dependence on the phase angle 0cr(k). 



A. First-Order Coupling 

The first-order subband coupling corresponds to direct mixing between \E'+(k) and ^E'-(k) 
and has straightforward description within nearly degenerate perturbation theory, where 
k = {kt, 0) for metallic SWNTs. Since 0o-(k) is continuous in vicinity of kf = 0, the diagonal 
matrix element M{kt,a;kt,(j) is continuous as well and merely shifts the location of the 
Fermi point and renormalizes the Fermi velocity. The change of the band gap is therefore 
determined by the off-diagonal term M{kt, a; kt, —a). 

According to Eq. ()30p. only Fourier components of the scalar potential with q = Q's 
contribute to the direct coupling. For example, for {n, n) armchair or (n, 0) zigzag nanotubes, 
it can be an angular perturbation of the form cos(2n^). Such Fourier component can be 
obtained by applying torsion, using chemical/biological decoration of the tube surface or the 
high multipoles of inhomogeneous potentials. On the other hand, since the pseudo-spinors 
of ^E'+(k) and ^l'-(k) are always orthogonal, the matrix element in Eq. ()30|) is reduced to 

Q ^ 



= ^[Ua- Ub) /2, (32) 

with Ua,b = Y1iqUq^^'^^°'^'°'^ ■ -^l- (E21) indicates that such potential components must be 
distinguishable at the two sublattices so as to mix the two orthogonal pseudo-spinors directly. 
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For the tensor potential, the symmetry is lower and even a uniform deformation (g = 0) 
can result in the first order coupling. The matrix element in Eq. ()31|) is reduced to: 

3 



M{kt,+;kt,-) = J] 57a cos arg(A:i) + - (k + K) ■ fa 



A=l 



kt^O 



sign(A;() [^71 cos^r] — 271/3) + 6-^2 cos r] + ^73 cos(r7 + 27r/3)] 

57i + 573\ , 



sign(/ci 



COST] I 572 — 



+ — sin ?7 (571 - 573) 



, (33) 



and more specifically for metallic achiral SWNTs, 



M(0,+;0, 



6'~f2 — ^73, armchair 
2(572 - (571 + 573), zigzag 



(34) 



which is consistent with previous findings about band gap changes in metallic SWNTs under 
uniaxial and torsional strain.— 



B. High-Order Coupling 

When the first order coupling between states ^^^-(k) and ^E'-(k) is forbidden, one has 
to turn to higher orders of the perturbation. Deriving the coupling could be tedious but 
some general rules can be built using the symmetry of pseudo-spinors. Similar to the case 
of A-SWNTs, coupling between the two crossing subbands of arbitrary metallic SWNT can 
be represented by H^^{kt) = ^^^{^j-y H^l{{'^i}), and the lowest contributing coupling 
order /iq can be determined by the dominating Fourier components of the potential. 

We first discuss the case of a scalar potential V{9) and restrict Q = in Eq. (jHUj) . The 
cancellation rule is similar to those which led to Eq.()16p. Assume there is a fi-th order 
coupling process between \l'+(k) and \I^-(k) through intermediate states = "^^^{^1,171^)} 
with i = 1 ■ ■ ■ n — 1: 



1=1 



cos 



{k, rfii^i) - {kt, rrii) 



i=l 



(35) 



where subscripts "0" and "/x" correspond to ^+(k) and ^^-(k) respectively, with mo = 
rh^ = 0. define a reversal process with intermediate states = ^-^^-^(^t, — m)}. After 
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rearrangement of the summation orders and using the relation (f)^cr{kt, ~rh) = 77 — (paikt, 1^), 
one arrives at 

H+^{{kt,-m^^i,-a^_i}) = — — 

^ {-ir-'H+4{kt,m,cTi}), (36) 

which cancel out with Eq. for even fi. The approximation sign for subband coupling pro- 
cess in Eq. (j36p arises from assumptions on the reflection symmetry E'^{kt, m) = E'^{kt, —rh) 
and reversal symmetry E'^{kt,rh) = —E^^{kt,rh), which can be derived from the linear 
dispersion approximation, i.e., E^{k) = crfF|k|. More generally, the reflection symme- 
try does not hold except for A-SWNTs due to the trigonal warping effect. For example, 
Ea-{kt, 1) 7^ E„{kt, —1) for metallic zigzag SWNTs, and a secondary band gap Eg oc al- 
ways opens under a uniform electric field perpendicular the nanotube radius.^ Here m = ±1 
is measured relative to the Fermi point K. When a nonzero overlap s or high order near- 
est neighbor interaction is included, the reversal symmetry can also be weakly broken. We 
conclude that the selection rules for arbitrary metallic SWNTs are similar to those for A- 
SWNTs, but may acquire a chirality dependence beyond the linear dispersion approximation. 

For general tensor potentials, M{kt;mi, ai;mj, aj) and M{kt, —rhj, —aj; —rhi, —ai) usu- 
ally have different magnitude and no simple cancellation rule can be built. An exception 
is when 72 = 73 for A-SWNT, e.g., under a radial deformation. In that case, the coupling 
between vr and n* subbands is reduced to zero as shown in Sec. IIVI due to cancellation from 
the reversal process within the nearest neighbor approximation. 

VII. CONCLUSION 

In this paper we employ the group theory approach to clarify the issue of MST in armchair 
and other metallic SWNTs under angular perturbations. We study symmetry requirements 
of MST on the nanotube and the perturbation potential, and demonstrate that the smallness 
of the MST effect is related to the symmetry of the pseudo-spinor components of the electron 
wave functions. Namely, the spinors of the crossing subbands are orthogonal, thus, any 
interaction between them is strongly weakened. For A-SWNTs, the gap is diminishing for 
almost any pure angular perturbation with a single angular Fourier component, due to the 
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XT ^^i—i—nii COS 




high couphng order. The couphng order is proportional to the number of atoms along the 
tube circumference for both types of perturbation studied here: the on-site (scalar) potential 
and off-site (tensor) deformation. The MST effect can be greatly enhanced by combining 
perturbations of different types and/or different angular momentums. 

Wc formulate selection rules for the band gap opening and its dependence on the per- 
turbation strength. The combination of the diagrammatic derivation of interaction matrix 
elements and group theory technique allows one to predict the scaling of the band gap on 
the potential: Eg oc |V|'^°, where the scaling exponent //q can be easily calculated for an 
arbitrary metallic SWNT for given symmetry of the potential. Corrections may arise due to 
refinement of the model TB Hamiltonian, e.g., the electron- hole asymmetry, the inclusion of 
(T-orbitals and others not considered here. As an example of such refined model, wc calculate 
the gap dependence on the overlap integral added to the classic orthogonal TB model. 

We present the analytical expression for the renormalization of the Fermi velocity, which 
occurs even if no MST is observed. The decrease of the Fermi velocity due to the perturbation 
is also seen as the enhancement of the DOS close to the Fermi level. 

The MST effect by the SWNT symmetry breaking could have potential applications for 
nanoscale electronic and optoelectronic devices. Additionally, we emphasize the possibility 
of engineering the nanotube DOS even when MST is forbidden under given perturbations, 
which can be potentially employed for SWNT opticals as well as switching devices. 
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APPENDIX: NEARLY DEGENERATE PERTURBATION THEORY 

When calculating the coupling between two nearly degenerate states and ^jj, it is 
more convenient to treat them as degenerate states. Since — Ej^ — only at the crossing 
point, one can include the energy dispersion by redefining the unperturbed Hamiltonian and 



18 



the external perturbation, for \Ec, — Ep\ being small: 

Ho = Ho- 1^,) - Ep 

Hi = Hi + Ho -Ho. (A.l) 

\1/q, and \&/3 now become degenerate states of Ho with £"0 ,3 = 0. Their original energy 
difference is absorbed in the perturbation while other states are not affected: 

^1 I ^,,;3) =£«,/?, (^a|^l|^,3) = (^a|i^l|^/3) 

(^'1 Hi 1^,,^) = (^'1 Hi 1^,,^) , (^'1 Hi = (^'1 Hi 1^") . 

As long as Hi — Hi is small, the rearrangement of Ho and Hi will not affect the results of 
the perturbation theory. 
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FIGURE CAPTIONS 



Fig. [It (color online). Schematics of second order coupling between n and n* subbands 
and the corresponding phase angles of intermediate states. 

Fig. [5J (color online). Band gap variation of (a) (5,5) and (b) (6,6) A-SWNTs as a 
function of the applied angular potential with q = 2. Insets: the unwrapped unit cell and 
schematics of the potential. Mirror planes of the potential pass through atomic sites so that 
all vertical mirror reflection and glide reflection symmetries are simultaneously broken. 

Fig. IHt (color online). Renormalized Fermi velocity vp of a (10, 10) A-SWNT as a function 
of u with q = 1 (circles) and q = 2 (squares). Solid and dashed lines are corresponding 
predictions from jQ(2u/q). Insets show the DOS structure near the Fermi level at m = 
and u = 1, with q = 1. 

Fig. 0t (color online). Band gap variations of A-SWNTs (a) as a function of u at e =0.05, 
0.1, and (b) as a function of e at m =0.5, 1. Dashed hues are to guide the eye. 
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FIG. 3: 
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